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Abstract 

Background: Proteins from the ABC family (ATP-binding cassette) represent the largest known group of efflux pumps, 
responsible for transporting specific molecules across lipid membranes in both prokaryotic and eukaryotic organisms. 
In arthropods they have been shown to play a role in insecticide defense/resistance. The presence of ABC transporters 
and their possible association with insecticide transport have not yet been investigated in the mosquito Anopheles 
stephensi, the major vector of human malaria in the Middle East and South Asian regions. Here we investigated the 
presence and role of ABCs in transport of permethrin insecticide in a susceptible strain of this mosquito species. 

Methods: To identify ABC transporter genes we obtained a transcriptome from untreated larvae of An, stephensi and 
then compared it with the annotated transcriptome of Anopheles gambiae. To analyse the association between ABC 
transporters and permethrin we conducted bioassays with permethrin alone and in combination with an ABC inhibitor, 
and then we investigated expression profiles of the identified genes in larvae exposed to permethrin. 

Results: Bioassays showed an increased mortality of mosquitoes when permethrin was used in combination with the 
ABC-transporter inhibitor. Genes for ABC transporters were detected in the transcriptome, and five were selected 
{AnstABCBl, AnstABCB3, AnstABCB4, Z\nsfABCmember6 and AnstABCG4). An increased expression in one of them 
(/AnsfABCG4) was observed in larvae exposed to the LD50 dose of permethrin. Contrary to what was found in 
other insect species, no up-regulation was observed in the AnstABCB genes. 

Conclusions: Our results show for the first time the involvement of ABC transporters in larval defense against 
permethrin in An. stephensi and, more in general, confirm the role of ABC transporters in insecticide defense. 
The differences observed with previous studies highlight the need of further research as, despite the growing 
number of studies on ABC transporters in insects, the heterogeneity of the results available at present does not 
allow us to infer general trends in ABC transporter-insecticide interactions. 
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Background 

Malaria is a major threat to human health and socio- 
economic development, representing a great burden in 
the vast regions of the world in which this parasitosis is 
endemic [1-3]. WHO estimated over 200 million cases 
of malaria in the 99 endemic countries and around 
660,000 deaths, in the year 2010 [2]. 
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Vector control through insecticides is a core compo- 
nent of malaria control programmes. However, continu- 
ous use of insecticides has led to the development of 
resistance in many malaria vectors around the world, 
which poses a serious threat to the global malaria con- 
trol efforts [3-5]. Research is therefore needed to under- 
stand the molecular basis of insecticide detoxification 
and develop even more effective methods to delay emer- 
gence of resistance [6]. 

In recent years, the role of ATP-binding cassette (ABC) 
transporters in the defense against toxic compounds as 
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pesticides has attracted a great deal of attention (reviewed 
in [7,8]). ABC transporters are ATP-dependent efflux 
pumps belonging to the ABC protein family located in the 
cellular membrane in both prokaryotic and eukaryotic or- 
ganisms. In eukaryotic organisms, they mediate the efflux 
of compounds from the cytoplasm to the outside of the 
cell or into organelles. ABC proteins have been subdivided 
into eight subfamilies (from ABC-A to ABC-H), and can 
transport a wide array of different substrates across cellu- 
lar membranes (e.g., amino-acids, sugars, lipids, and pep- 
tides). Most of the ABC transporters associated with the 
efflux of pesticides belong to the subfamilies ABC-B (also 
referred to as P-glycoproteins, P-gps), ABC-C and ABC- 
G. In some cases ABC-transporter action has also been 
associated with insecticide-resistant phenotypes in species 
of agricultural or medical importance [7,8]. In spite of an 
increasing awareness of the potential importance of ABC 
transporters in vector control, to date they have been 
poorly studied in detail in malaria vectors [8] . Here, we in- 
vestigated the role of ABC transporters in the detoxifica- 
tion against the insecticide permethrin in the malaria 
vector Anopheles stephensi (Culicidae: Diptera). This 
mosquito species, vector of both Plasmodium falcip- 
arum and Plasmodium vivax, is one of the major vec- 
tors of human malaria in the world. An. stephensi 
occupies a geographic range that spans from the Mid- 
dle East to South-East Asia [9]. These regions contrib- 
ute to 15% of malaria cases worldwide, with an estimated 
28 million people annually affected by the disease [2]. Per- 
methrin belongs to the pyrethroid class of insecticides, 
which is by far the most commonly used in malaria 
vector-control interventions [2]. 

Pyrethroids act by modifying the gating kinetics of 
voltage-gated sodium channels, thereby disrupting neuron 
function, which leads to rapid paralysis and death of the 
insect [10]. They can enter into the insect body by inges- 
tion and penetration into the hemolymph through the ali- 
mentary canal, or via contact with sensory organs of the 
peripheral nervous system [11]. Insect midgut is rich in 
ABC transporters, whose action, therefore, likely prevents 
permethrin to reach its target sites [8]. Furthermore, in- 
sects possess protective neural barriers (e.g. a layer of gli- 
ally derived epithelial cells), where ABC transporters likely 
play an important role in the exchange of molecules 
[12,13]. In particular, inhibition of P-gp in Schistocerca 
gregaria has been shown to increase brain uptake of differ- 
ent drugs [13]. The involvement of ABC transporters in 
pyrethroid detoxification has been reported for a few in- 
sect species, such as Helicoverpa armigera [14-16], Apis 
mellifera [17] and Culex pipiens [18]. Up-regulation of 
ABC-transporter genes has also been reported in pyrethroid 
resistant strains of the bed bugs Cimex lectularius [19] and 
of the vector mosquitoes Anopheles gambiae [20] and Aedes 
aegypti [21]. No protein belonging to the ABC transporters 



has yet been described in larvae of An. stephensi, nor the 
possible association of this class of proteins with insecticide 
transport has been investigated in this species. In this paper 
we investigated the presence and role of ABCs in transport 
of permethrin insecticide in larvae of An. stephensi: i) by bio- 
assays with permethrin alone and in combination with an 
ABC inhibitor; ii) by investigating gene expression profiles 
in larvae exposed to permethrin treatment. 

Methods 

Mosquito samples 

The mosquito larvae used in this study were obtained 
from adult females of a An. stephensi colony, derived 
from the Liston strain. This colony has been maintained 
for four years in the insectary at the University of 
Camerino, following standard conditions: adult insects 
are reared at 28 ± 1°C and 85-90% relative humidity with 
photoperiods (12:12 L-D) with a 5% sucrose solution, 
and adult females are fed with mouse blood for egg lay- 
ing. Eggs from this colony were put into spring water in 
order to obtain the larvae. Larvae were maintained in 
spring water and fed daily with fish food (Tetra, Melle, 
Germany) under the same conditions as the adults. 

Bioassays 

Inhibition of ABC-transporters should lead to a higher 
intracellular concentration of insecticide, thus increasing 
larval susceptibility and insecticide efficacy [8]. In order 
to evaluate a potential synergy, we performed bioassays 
with permethrin insecticide alone and with permethrin 
in combination with a sub-lethal dose of the ABC- 
transporter inhibitor verapamil (see below for experi- 
mental determination of sub-lethal dose of verapamil). 
This is a calcium channel blocker, which works by com- 
peting with cytotoxic compounds for efflux by the mem- 
brane pumps [22]. All bioassays were conducted on An. 
stephensi larvae at the third instar, according to standard 
protocols [23]. 

Groups of 25 larvae were put in 250 ml plastic glasses 
with 100 ml of spring water and different concentrations 
of insecticide or insecticide + inhibitor. All tests were 
performed in quadruple. Additional groups of larvae, 
treated only with water and acetone (that was used to di- 
lute permethrin), were used as controls. Mortality was 
assessed at 24 h post-treatment and the larvae were con- 
sidered dead if immobile, even after a mechanical stimulus. 

In the bioassays with permethrin alone (Sigma- Aldrich 
S.r.l., Milan, Italy), six insecticide concentrations were 
used (0.015, 0.047, 0.092, 0.23, 0.57, 1.44 mg/1) to have 
mortality in the range 1-99%. The drug was dissolved in 
acetone and then diluted in water to obtain the test so- 
lutions. The bioassays with permethrin in combination 
with verapamil were performed using permethrin at the 
six concentrations indicated above, plus two additional 
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concentrations (0.0024 and 0.0048 mg/1). The sub-lethal 
dose of verapamil (i.e. the dose at which no dead larvae 
were observed) was determined using ten different con- 
centrations (20, 40, 80, 100, 160, 240, 320, 400, 480, 
560 uM) following the protocol above. The larval mor- 
tality data were subjected to Probit regression analysis 
[24] as implemented in the XLSTAT-Dose software (avail- 
able at: http://www.xlstat.com) to estimate the LD50 
values and their 95% confidence intervals (CIs). To es- 
timate the effect on larval mortality of the ABC inhibi- 
tor at sub-lethal dose, the synergistic factor (SF) was 
calculated. 

Identification of ABC transporter genes 

A total of 200 untreated larvae of An. stephensi at the third 
instar were pooled in 15 ml of RNAlater stabilization solu- 
tion (Qiagen, Hilden, Germany) and provided to an exter- 
nal company (GATC Biotech AG, Costance, Germany) for 
one run of 2x250 paired-ends reads sequencing on the Illu- 
mina MiSeq platform. The resulting reads were assembled 
using Trinity with default settings [25]. The assembled 
contigs were compared with Blastx (evalue 0.00001) to the 
annotated transcriptome of An. gambiae available in the 
VectorBase database, and the sequences of ABC trans- 
porters were extracted automatically and manually con- 
trolled. Based on published results about the involvement 
of ABCs on multidrug resistance in several arthropods 
(mainly mosquitoes) [8], we selected five genes from the 
transcriptome of An. stephensi. Oligonucleotide primers 
were then designed from the sequence of each gene 
(Table 1). The sequences of ABC transporters identified in 
An. stephensi were translated to aminoacids and compared 
against the UniProt database [26] using Blastp. Homolo- 
gous proteins were aligned using ClustalX [27] and dis- 
tances among them were estimated by Dayhoff PAM 
matrix as implemented in the PROTDIST software of the 
PHYLIP package [28]. 

Gene expression profile after insecticide treatment 

The activity of ABC-transporters is generally modulated 
at gene transcriptional level: the presence of toxic com- 
pounds leads to higher transcription. In order to assess 
this topic, larvae of An. stephensi at the third instar were 



exposed to permethrin and the expression of ABC- 
transporter genes was monitored in the surviving larvae 
by quantitative RT-PCR twice after insecticide treatment: 
24 h (e.g. the time at which the LD50 has been esti- 
mated) and 48 h following the study of Figueira-Mansur 
[29] that found increased expression of ABC trans- 
porters in the mosquito Ae. aegypti. The larvae were 
treated with the LD50 (0.137 mg/1) of insecticide esti- 
mated by bioassays as described above and two pools, of 
ten larvae each, were collected after 24 and 48 h of 
insecticide-treatment. All pools of larvae were stored in 
RNAlater for molecular analysis and, controls (water + 
acetone) were collected following the same time frame. 

RNA was extracted from each pool of larvae using the 
RNeasy Mini Kit (Qiagen, Hilden, Germany) including an 
on-column DNase I treatment (Qiagen, Hilden, Germany), 
according to the manufacturer's instructions. Total RNA 
was eluted into nuclease-free water and the concentration 
of RNA was determined at 260 nm [30] using a NanoDrop 
ND-1000 (Thermo Scientific, Delaware, USA). cDNAs 
were synthesized from 250 ng of total RNA using a 
QuantiTect Reverse Transcription Kit (Qiagen, Hilden, 
Germany) with random hexamers. The cDNA was used 
as template in RT-PCR reactions using the primers de- 
signed from the sequences of identified ABC genes 
(Table 1). The amplification fragments, obtained using 
standard PCR conditions and the thermal profile indi- 
cated below, were sequenced in order to confirm the 
specificity of the amplification. 

Quantitative RT-PCRs on the target ABCs were per- 
formed using a BioRad iQ5 Real-Time PCR Detection 
System (Bio-Rad, California, USA), under the following 
conditions: 50 ng cDNA; 300 nM of forward and reverse 
primers; 98°C for 30 sec, 40 cycles of 98°C for 15 sec, 
59"C for 30 sec; fluorescence acquisition at the end of 
each cycle; melting curve analysis after the last cycle. 
The cycle threshold (Ct) values were determined for 
each gene, in order to calculate gene expression levels of 
target genes relative to rps7, the internal reference gene 
for An. stephensi [31]. The expression of the ABC trans- 
porters genes in the control group was considered as the 
basal level (equal to 1). The estimates of the expression 
level of each gene in the treated larvae are reported as 



Table 1 Primer sequences of ABC transporter genes identified in Anopheles stephensi 


Gene 


Forward primer 


Reverse primer 


PCR product size (base pairs) 


Source 


/\n5fABCB2 


TATCAAGTTCACGGATGTAGAGT 


TATCCACCTTGCCACTGTC 


185 


This work 


/\n5fABCB3 


CAACCGTTCCGTAATACTACC 


ACTGGTAGCCCAATGTGAAG 


133 


This work 


AnstN&C&A 


GGACAAAACATTCGGGAGG 


CGTAGTGAATGTTGTGGCG 


109 


This work 


/\n5fABCBmemb6 


CTGGAGACGCTGAGAGATA 


TACTCCTCGGTGAACTGG 


125 


This work 


/\ri5fABCG4 


ATGAGCCCATTCGTCCTG 


AGCGTGGAGAAGAAGCAG 


158 


This work 


Rps7 


AGCAGCAGCAGCACTTGATTTG 


TAAACGGCTTTCTGCGTCACCC 


90 


Capone et al. 2013 [31] 
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the means ± standard deviation (SD) in Additional file 1: 
Table SI. 

Ethical statement 

Maintenance of the mosquito colony of An. stephensi 
was carried out according the Italian Directive 116 of 
10/27/92 on the "use and protection of laboratory ani- 
mals" and in adherence with the European regulation 
(86/609) of 11/24/86 (licence no. 125/94A, issued by the 
Italian Ministry of Health). 

Results 

Bioassays 

No mortality was observed when larvae were exposed at 
concentrations of verapamil up to 100 uM; this concen- 
tration was thus used as the sub-lethal dose in the bioas- 
says with insecticide + ABC-transporter inhibitor. The 
results of toxicity assays using permethrin and permeth- 
rin in combination with verapamil are reported in 
Table 2. The mortality data observed in bioassays well 
fitted the Probit dose-response model (Chi-Square prob- 
ability <0.0001). The LD50 dose in permethrin assay was 
0.137 mg/l while in the assay in combination with verap- 
amil LD50 was 0.025 mg/l (Table 2). No overlapping 
values were observed between LD50 95% CI of insecti- 
cide alone and insecticide plus verapamil; the addition of 
verapamil increased the toxicity of permethrin about 
5-fold (SF = 5.48). 

Isolation of ABC transporter genes and expression profile 
after insecticide treatment 

The Illumina MiSeq platform was used to sequence the 
cDNA library obtained from a pool of 200 larvae of An. 
stephensi, and 16,686,276 paired-ends reads were ob- 
tained. MiSeq raw data were assembled with Trinity, 
obtaining 40,498 contigs. The contigs containing ABC 
transporter genes were extracted on the basis of the an- 
notated transcriptome of An. gambiae available in data- 
base. Five sequences, respectively of 3612, 2154, 2481, 
2553 and 2182 base pairs, were found to share 85-94% 
identity with putative ABC multidrug transporters of An. 

Table 2 Toxicity of verapamil and permethrin against 
Anopheles stephensi larvae 

Insecticide N Slope ± SE LD50 (95% CI) SF 

Verapamil 600 3.846 (0.374)* 528 \iM (486-587) 

Permethrin 600 1.819 (0.125)* 0.137 mg/l 

(0.117-0.160) 

Permethrin + verapamil 600 2.123(0.174)* 0.025 mg/l 5.48 
(100 uM) (0.021-0.029) 

LD50 and slopes of regression lines estimated from mortality data by Probit 
analysis are shown. N, number of larvae used in bioassays; SE, standard error; 
95% CI, 95% confidence interval. SF, synergistic factor. 
•Chi-Square probability < 0.0001. 



gambiae: ABCB2 (AGAP005639) (85% identity), ABCB3 
(AGAP006273) (94% identity), ABCB4 (AGAP006364) 
(88% identity), ABCmember6 (AGAP002278) (94% iden- 
tity) and ABCG4 (AGAP001333) (85% identity). We de- 
noted them as AnstARCKl, A«stABCB3, AnstABCM, 
A«sfABCmember6, ^4«s£ABCG4 and we deposited them in 
EMBL Nucleotide Sequence Database [EMBL: LK392613 
to LK392617]. The alignment of the deduced amino acidic 
sequences of ABC transporters identified in An. stephensi 
with sequences of homologous ABC transporters of An. 
gambiae is shown in Additional file 2: Figure SI. Dayhoff 
PAM distance estimates between the ABC transporters 
identified in An. stephensi and the homologous ABC trans- 
porters of mosquitoes An. gambiae, Anopheles darlingi, 
Ae. aegypti and Culex quinquefasciatus that showed the 
highest percentage of identity following Blast search are 
shown in Additional file 3: Table S2. 

Conventional PCR amplicons obtained from each gene 
primer set were sequenced, confirming in all cases the 
sequences generated with the MiSEQ experiment The 
RT-PCRs were performed to investigate whether per- 
methrin treatment at the LD50 dose (0.137 mg/l, Table 2) 
increased or decreased the ABC gene expression in the 
An. stephensi larvae after 24 and 48 h of insecticide- 
treatment. As reported in Figure 1 and Additional file 1: 
Table SI, after 24 h of permethrin treatment, the relative 
expression of all selected genes was down-regulated, 
with the exception of the AnstABCG^ gene, that showed 
about three-fold increase of expression compared to the 
control. Similarly, after 48 h of permethrin treatment, 
the relative expression of all ABCB genes was down- 
regulated, while the A«siABCG4 gene showed a ten-fold 
increase of expression compared to the control. 

Discussion 

Bioassays and molecular data suggest the involvement of 
ABC transporters in the defense of An. stephensi larvae 
against the permethrin insecticide. Indeed, inhibition of 
ABC-transporters led to a higher susceptibility of larvae 
to insecticide, indicating that ABC transporters are asso- 
ciated with insecticide detoxification. In addition, in 
mosquito larvae exposed to the LD50 dose of permeth- 
rin, we observed an increased expression of AnstNRCGA, 
one of the five tested genes coding for ABC transporters. 

Arthropod ABCG genes are orthologous of the human 
gene ABCG2, which has been associated with resistance 
to cancer drugs [32,33], while data in insects show that 
ABCG transporter genes were significandy over-transcribed 
in response to exposure to insecticides. Microarray gene ex- 
pression studies revealed that ABCG transporter genes were 
up-regulated in DDT resistant strains of Drosophila melano- 
gaster [34] and in a Plutella xylostella (Lepidoptera) strain 
resistant to chlorpyrifos [35]. The ABCG4 transporter gene 
was over-transcribed in Bemisia tabaci whiteflies resistant to 



Epis et al. Parasites & Vectors 2014, 7:349 
http://www.parasitesandvectors.eom/content/7/1/349 



Page 5 of 7 




0.01 



100 



10 



48 h 



0.1 



0.01 



ABCB2 ABCB4 ABCG4 

ABCB3 ABCB6 

Figure 1 Relative expression of Anopheles stephensi ABC genes 
measured by quantitative PCR after 24 and 48 h of permethrin 
exposure. The expression level in non-treated larvae was considered 
to be the basal level (equal 1). The internal reference gene rps7 for 
An. stephensi was used to normalize expression levels. The values are 
expressed as means ± standard deviations. 



the neonicotinoid thiamethoxam [36] as well as in Anoph- 
eles arabiensis resistant- and sensible-strains to DDT [33]. 
The ABCG3 gene was found differentially expressed in Ae. 
aegypti pyrethroid resistant populations versus susceptible 
strains [21]. 

The other four genes coding for ABC transporters that 
we detected and tested in An. stephensi (AnstABCB2, 
AnstABCBZ, A«sfABCB4 and AKs£ABCmember6) be- 
long to the ABCB subfamily. Several members of this 
subfamily have been associated with transport and/or re- 
sistance to different insecticide classes in several insect 
species [7,8]. In mosquitoes, the ABCB2 gene was 
showed by quantitative PCR to be eightfold up-regulated 
in larvae of a susceptible strain of Ae. aegypti analysed at 
48 h post temephos-treatment [29]. The ABCB4 gene was 
showed to be over-transcribed by transcriptome and 
quantitative PCR analyses, in DDT and pyrethroid (per- 
methrin and deltamethrin) resistant Ae. aegypti popula- 
tions compared to a laboratory susceptible population 
[21]. Microarray analysis showed that the ABCB4 was up- 
regulated in different populations of DDT-resistant An. 
gambiae mosquitoes [37]. Our results showed no over- 



transcription of the ABCB genes in An. stephensi suscep- 
tible larvae exposed to permethrin, while insecticide treat- 
ment induced an increased expression of the ^4«s£A.BCG4 
(Figure 1, Additional file 1: Table SI). 

On the whole, the results herein presented support the 
view of the involvement of ABC transporters in insecti- 
cide transport, although differences with previous stud- 
ies have been observed. Are these differences due to the 
insecticides used or to the status of the analysed samples 
(i.e. susceptible vs. resistant)? Despite the growing num- 
ber of studies on ABC transporters in insects, the het- 
erogeneity of the data available at present does not allow 
to infer general trends that may underlie particular inter- 
actions between ABC transporters and insecticides. Fur- 
ther studies are needed to highlight these and other 
issues. For example, our results showed that the expres- 
sion of ABCB genes in An. stephensi did not only in- 
crease in larvae treated with permethrin compared to 
those non-treated, but indeed it decreased (Figure 1, 
Additional file 1: Table SI). In the latter case, the study 
of gene expression at more time points could contribute 
to the understanding of whether there are temporal de- 
lays, or whether compound-specific or species-specific 
differences exist in their activation [38-40]. Furthermore, 
most studies have been conducted on larval stages [7,8]. 
The synergist and transcript profiles may differ between 
larval and adult stages, an interesting topic for future 
studies. 

Diffusion of vectors of human diseases driven by hu- 
man activities and global climate change as well as in- 
surgence of insecticide resistance can seriously impact 
our ability to control vector-borne diseases [41-44]. Fur- 
thermore, environmental pollution and resistance phe- 
nomena clearly show the limits of the chemical approach 
for pest control and the need to delineate new strategies 
that optimize the use of available molecules, with the aim 
of reducing their impact on the environment [31,45-48]. 

In the last decades advances in molecular techniques 
have greatly improved our tools to investigate the dy- 
namics of vector populations and of pesticide resistance 
insurgence [43,49-53]. More recently, next-generation 
sequencing technologies have offered unprecedented op- 
portunities to investigate the molecular basis of the 
interaction between cellular defenses and insecticides 
[54]. In this context, the increasing interest about ABC 
transporters in transport and/or resistance against insec- 
ticides led to an increase of the information on these 
genes in various insect species and their association with 
insecticide detoxification [7,8]. 

Conclusions 

In this study we have demonstrated for the first time in 
the larvae of An. stephensi that verapamil increases the 
sensitivity to permethrin in laboratory assays; in addition, 
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we isolated five genes encoding for ABC transporters, and 
investigated their expression profile after exposure to per- 
methrin. To analyse the potential role of ABC transporters 
in permethrin transport in An. stephensi, we performed 
bioassays using a sub-lethal dose of the ABC transporter 
inhibitor verapamil in association with permethrin. The 
results obtained using this approach highlight that the 
combination of insecticides with an ABC-transporter in- 
hibitor can increase the efficacy of the insecticide mol- 
ecule [18,29,55]. In prospect, combined treatments of 
insecticide plus ABC-transporter inhibitors could be pro- 
posed, with the objective of reducing the current dosages 
of insecticides or to prevent the development of resist- 
ance, and reduce environmental pollution [29,56]. The im- 
plementation of such a strategy would require the availability 
of gene- and species-specific inhibitors in order to avoid the 
serious consequences that would derive from a generic in- 
hibition of ABC-transporters in non-target organisms. The 
study of ABC-transporters at the gene level is therefore cru- 
cial for the understanding of both their potential role as 
defense mechanisms and for their inhibition for vector con- 
trol purposes. 
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